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Following the simple observation that the interconnection of a set of quantum op- 
tical input-output devices can be specified using structural mode VHSIC Hard- 
ware Description Language (VHDL), we demonstrate a computer-aided schematic 
capture workflow for modeling and simulating multi-component photonic circuits. 
We describe an algorithm for parsing circuit descriptions to derive quantum equa- 
tions of motion, illustrate our approach using simple examples based on linear and 
cavity-nonlinear optical components, and demonstrate a computational approach 
to hierarchical model reduction. 
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Ongoing advances in materials science, nanoscale synthesis and lithography are es- 
tablishing key technical requirements for the fabrication of complex nanophotonic 
circuits. Critical applications for such circuitry can be foreseen in ultra-low power 
photonic signal processing and interconnects, sensor networks, and quantum in- 
formation technology. Although current research in nanophotonics focuses largely 
on the physics of individual devices such as resonators and waveguides, there are 
substantial new challenges to be addressed in the analysis and design of circuits 
comprising large-scale networks of interconnected components with dynamic opti- 
cal nonlinearities. For premier applications, quantum-optical network models will 
be required in order to capture fluctuations, coherence, and entanglement effects 
that may be decisive drivers of the overall circuit performance. Suitable theoretical 
frameworks exist but require cumbersome algebraic manipulations for the derivation 
of multi-component models, pointing to the need for computer-aided paradigms for 
generating quantum-optical equations of motion from high-level circuit representa- 
tions that can be manipulated more intuitively by circuit engineers. Simply put, we 
arc rapidly approaching a new phase of research in quantum nonlinear photonics 
in which we will need user-friendly circuit design tools like the ones we have long 
exploited in classical electronics. 

In this paper we propose and demonstrate a modeling and simulation work- 
flow based on schematic capture using a Quantum Hardware Description Language 
(QHDL) for nanophotonic circuits, which we will define as a proper subset of the 
standard VHSIC Hardware Description Language (VHDL). Our approach utilizes 
a mixture of common open-source software packages and custom processing scripts 
to provide a high-level, modular interface to the quantum circuit 'algebra' of Gough 
and James [l|, 0] (which generalizes earlier work on cascaded open quantum systems 
by Carmichael Q and by Gardiner Q). The natural hierarchical organization of 
VHDL and the schematic capture workflow should facilitate future work on model 
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reduction and design abstractions for nanophotonic circuits, which seems essen- 
tial given the extremely high dimension (variable count) associated with many- 
component quantum models. 

In the following sections we first review the formal setting of (S, L, iJ) compo- 
nent models and the concatenation and series and products as introduced by Gough 
and James, which have recently been used to derive quantum nonlinearphotonic 
circuit models by hand or using custom-coded computer algebra scripts [aSH,!!!. 
While we will restrict our attention here to linear and cavity nonlinear optics, 
it should be noted that the (S,L,iJ) formalism can in principle be used to de- 
scribe hybrid circuits incorporating suitable spintronic, nanomechanical and/or 
quantum-electronic components. Likewise, the approach we describe here could be 
extended straightforwardly to admit static Bogoliubov components as described 
in [l3] ■ We then review the proposed syntax of QHDL and illustrate its use in the 
specification of a simple interferometer as a network of elementary optical compo- 
nents. After describing methods that can be used to parse QHDL circuit descriptions 
to derive quantum equations of motion for analysis and/or numerical simulation, 
we illustrate the full schematic capture workflow using an example of constructing 
a bistable latch from cavity nonlinear optical components. The paper closes with a 
brief consideration of model reduction in the (S, L, H) context. 



1. Modeling quantum circuitry 

Within this section we use {Qj, j ~ 1,2,3,..., N} to denote individual quantum 
input-output components. We clearly distinguish between input and output ports 
and do not consider bi-directional ports, although for physical reasons every input 
port is assumed to have an associated output port and vice versa. 



(a) The circuit algebra 

Our modeling workflow is based on the Gough- James synthesis results for open 
quantum systems [l|, 0] , which provide a purely algebraic method to derive quantum 
Markov models for a network of interconnected quantum components. 

A component with an equal number n of input and output channels is described 
by the parameter triplet (S, L, H), where H is the effective internal Hamilton oper- 
ator for the system, L = (Li,i2, ■ • • ,L„)'^ the coupling vector and S = [Sjk)^^^-^ 
is the scattering matrix^ whose elements are themselves operators. 

Each element Lk of the coupling vector is given by an operator that describes 
the system's coupling to the fc-th input channel. Similarly, the elements Sjk of the 
scattering matrix are given by system operators describing the scattering between 
different field channels j and k. The only mathematical conditions on the parameters 
are that the Hamiltonian operator be self-adjoint and the scattering matrix be 
unitary: 

iJ^ = iJ and S^S = SS^ = 1„. 
The master equation [ll| corresponding to a given (S, L, H) model is 

^ = -^[H,p,] + (l,p,L* \ {L*L,,p,]) (1.1) 
i=i ^ ^ 
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Here [A, B] = AB - BA and {A, B} = AB + BA, while pt is a density matrix 
describing the evolving state of the internal degrees of freedom. It is also straight- 
forward to obtain the quantum filtering equations [l^, [I^ for stochastic simulation 
of a given (S, L, H) model. 

While the scattering matrix elements Sjk do not appear in Eq. (|l.ip they are 
required for the composition rules described below, which can be used to derive the 
overall parameter triplet for a network of interconnected quantum input-output 
components. The (S, L, H) circuit algebra plus simple correspondences such as 
Eq. provide all that is needed to obtain overall equations of motion for complex 
photonic circuits. 




(a) Qi ffi Q2 (b) Q2 < Qi (c) [Q]k^i 

Figure 1. Basic operations of the Gough- James circuit algebra. 



In [2[, Gough and James have introduced two operations that allow for the 
construction of arbitrary networks of optical feedforward circuits: 



1. The concatenation product (cf. Figure 1(a) I describes a formal adjoining of 



two systems in which there is no optical scattering between the systems: 



(Si, Li, i/i) ffl (S2, L2, -ff2) 



Si 
S2 



(1.2) 



Note however, that even without optical scattering, the two subsystems may 
interact via shared quantum degrees of freedom. A simple example of this 
scenario is given by a two-mode resonator (such as a ring-resonator) with an 
atom that interacts with both optical modes, but in which there is no direct 
scattering between the modes. 



2. The series product (cf. Figure 1(b) I describes a configuration in which two 
systems Qj = (Sj,Lj,i/j), j = 1,2 possessing an equal number of channels 
n are connected in such a manner that all output channels of Qi are fed into 
the corresponding input channels of Q2- The resulting system is then given 

by 

(S2, L2, iJ2) < (Si, Li, Fi) = (S2S1, L2 + S2L1, Fi + F2 + 3 {4S2L1 }) , 

(1.3) 



where we define the imaginary part of an operator as 



— 2i ■ 



To make the network operations complete, one additional rule is required: The 
feedback operation (cf. Figure 1 (c) I describes the case where the fc-th output channel 
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of a system with n > 2 channels is fed back into the ^-th input channel. The result 
is a system with n — 1 channels: 

[(S,L,iJ)],^,^ (s,L,ff) (1.4) 

Formulae for the resulting parameter triplet are provided in [Appendix A[ 

Note that the series product can be expressed in terms of the concatenation 
and feedback operations {e.g., for two components with n = 1 we have Q2 <i Qi = 
[Qi ffl(52]i-i.2), and consequently, the latter two operations arc sufficient to perform 
all network calculations. However, the series product is a useful shorthand and 
allows for a more intuitive network expression. 

For use in the following we define the identity system with n channels 

1„ = (1„,0,0), (1.5) 

where 1„ ~ ('5fci)I- 1=1 the identity matrix in n dimensions, as well as the channel 
permuting system 

P,, = (P,,0,0), (1.6) 

where the permutation matrix is defined by Po- = ('5fc.cr(i))^ ^ .|^. This definition 
ensures that Pg-, < Pen = Pa^ocn ■ 



(6) The QHDL syntax 



QHDL is a subset of structural VHDL [lj|, which we will use as a formal syn- 
tax for specifying photonic circuits in terms of interconnections among referenced 
quantum input-output components. These components can themselves represent 
composite networks of subcomponents, facilitating hierarchical approaches to pho- 
tonic circuit design. It is useful to start with a set of basic components such as 
beamsplitters and phase-shifts, as well as linear and non-linear cavity models with 
one or more coupling mirror^, which can be collected in a shared library file. The 
set of such primitive components within a QHDL software environment can of 
course be extended at any time. 

Within the context of a single QHDL file, the exact physical model (parameter 
triplet) of any referenced component is left unspecified except for its external ports 
and parametric dependencies. This approach allows the circuit designer to operate 
at a high level of abstraction, facilitating last-minute substitution of alternative 
physical component models (including effective models with reduced simulation 
complexity) into a given interconnection topology. 

In the following section, we will introduce the QHDL syntax by means of a very 
simple circuit that realizes a Mach-Zehnder interferometer. 

A QHDL file begins with the entity declaration, which defines the abstract 
interface of the circuit being specified: it specifies a list of named input and output 

t It is important to note tiiat the Gough-James circuit algebra cannot be used to build dynam- 
ical systems from static components, e.g., it cannot create the Fock space and operator algebra 
for an optical resonator mode as an automatic result of cascading beamsplitters and phase-shifts 
in the configuration of a ring cavity. All such dynamic components therefore must be implemented 
as primitive (S, L, H) models. 
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Figure 2. A basic Mach-Zehnder setup. 



ports (of the overall circuit), which arc required in order for the circuit itself to 
be callable as a composite QHDL component, as well as any numeric parameters 
required for physical modeling. Note that we require that all input ports appear 
before all output ports. 

Listing 1. Entity declaration 

entity Mach_Zehnder is 

generic (phi_mz: real := 0); 

port (Inl, Vacin: in ficldmodc; Outl, Out2: out ficldmode); 
end Mach_Zohndor; 

For this entity we must then have one or more architecture declarations in the 
same QHDL file. These provide alternative ways of realizing the internal structure 
of the circuit. The architecture declaration consists of a head which specifies the 
interfaces of all components used in the architecture body and all internal signals. 
The component declarations are very similar to the entity declaration- they serve 
to establish an interface for each subcomponent. 

Listing 2. Architecture head 

architecture structure_MZ of Mach_Zchnder is 
component beamsplitter 

port (a, b: in fieldmode; c, d: out fieldmode); 
end component beamsplitter; 

component phase 

generic (phi: real); 

port (a: in fieldmode; b: out fieldmode); 
end component phase; 

signal bsl_phase, bsl_bs2, phase_bs2: fieldmode; 

The architecture body then consists of a series of instance assignments for each 
occurrence of any of the previously specified component types. Each instance assign- 
ment specifies the relationship between the component-instance parameters and the 
entity parameters. In addition, it specifies a port map detailing how the component- 
instance is connected to the internal signals or the external ports. 

Listing 3. Architecture body 

begin 

BSl: beamsplitter 

port map (a => Inl, b => Vachi, c => bsl_bs2, d => bsl_phase); 
phase: phase 

generic map (phi => phLmz); 
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port map (a => bsl_phasc, b => phasc_bs2); 
BS2: beamsplitter 

port map (a => phaso_bs2, b => bsl_bs2, c => Outl, d => Out2); 
end structurc_MZ; 

In the port map, each internal component port is assigned to either an entity 
port or a signal. Any instance in (out) port must be connected either to an entity 
in (out) port or to a signal that is connected to another instance's out (in) port. 

Listing 4. Port map statement 
port map (a => phase_bs2, b => bsl_bs2, c => Outl, d => Out2); 

Each signal therefore connects exactly two ports: one instance input and one in- 
stance output or one instance input (output) and an entity input (output). 

(c) Parsing a network 

Here we present a simple algorithm to parse a general network into a circuit 
expression. We assume that the QHDL file has been preprocessed such that we 
have the lists of ports, components, instances, signals and port mappings in native 
data structures accessible to our algorithm. 

1. We denote the list of internal signals by 5*. For each instance assignment 
j = 1, 2 . . . iV in the architecture body: 

• Generate the network triplet Qj ~ (Sj , Lj , Hj ) with the correct parametriza- 
tion as specified in the generic map statement. 

• Generate the correctly ordere^l list of input port names Ij and the 
correctly ordered list of output port names Oj where each portname is 
entry is of the form instance-name:port-name. 

2. Concatenate all triplets Q ^ QiSiQ2Si- ■ -SiQ^ and similarly concatenate the 
input and output port lists I = Ii + 12 + ■ ■ ■ + In and O = Oi +O2 + ■ ■ ■ +On 

3. For each internal signal s € S* concatenate the full circuit triplet Q with a 
single channel identity system li resulting in Q^p — Q ffl 1|5|, 

4. Now, each element in the full list of output ports O corresponds to an entry 
of the form instance-name:port-name. Make copies of O' = O and S' ~ S and 
iterate over all output ports in the following fashion: 

If the output port is connected to a global output (i.e. an entity output port) , 
continue to the next entry. 

If the output port is connected to the j-th signal in the current signal list S", 
let k be the index of the output port in the current output port O' list and 
update the model triplet Q^"^ — >• Qy""*"^' = [Q^p^]k^M+j, where M — \0'\ is 
the length of the current output port list. Then, remove the A:-th entry from 
O' , and the j-th entry of S' . 

t As defined via the the component declaration in the architecture head. 
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5. Now, let Mf = \S'\ and iterate over a copy of the input port list I' = I and 
a new copy of the signal list S" = S: 

If the input port is connected to a global input (i.e. an entity input port), 
continue to the next entry. 

If the input port is connected to the j-th entry of S"' update Qj"^ Q^P'^^^ = 

[Q/"'']Mj:+j-i.fc where k is the index of the current port in /'. Then, remove 
the fc-th entry of /' and the j-th entry of S" . 

6. By construction, the only remaining ports of our resulting triplet lead 
to global/entity ports. Iterating over O' and the list of entity output ports 
Oe, construct a suitable permutation (Tout that maps every output port index 
from O' to the correct index of the entity output port. In a similar fashion, 
iterate over /' and the list of entity input ports Ie to generate a permutation 
(Tj~^, mapping the indices from /' to the correct indices of the entity input 
ports within Ie- Then, invert this permutation ain = (^in^)"^ to obtain a 
mapping from Ie to /'. Finally, the model triplet for the circuit is given by 

If one is interested in working with the actual network expressions as opposed 
to the more concrete level of the actual Hilbert space operators, there exist other, 
more complex, approaches to parsing a network, which directly yield simpler over- 
all network expression. Combined with a sufficiently sophisticated set of circuit 
expression simplification rules, the above algorithm works just as well. 

(d) The QHDL workflow 

The circuit design workflow relies heavily on symbolic computer algebra meth- 
ods. Using symbolic algebra, rather than working with numerical matrix represen- 
tations of all the operators appearing in the component parameter triplets, makes 
it possible to view the overall circuit (S,L,iJ) in analytic form. It also allows the 
designer to defer choosing the values of numerical parameters, which could be conve- 
nient for optimization scenarios, as well as details such as the upper photon-number 
limits to use for truncated Fock spaces in numerical simulations. 

In fact we can define our own algebraic types, operations and simplification 
rules not just for Hilbert space operators and scalar coefficients, but also for circuit 
algebra components. This approach enables us to extend the hierarchical design 
principle even to our compiled QHDL component library, as will become clear in 
the following outline of the modeling workflow: 

1. Circuit design In step 1, we visually compose the circuit using a schematic 
capture tool and then export to QHDlHI or directly describe the circuit in 
text-based QHDL. Since QHDL describes the connections between functional 
entities, it is not necessary at this stage to specify how referenced components 
are implemented. 

t In our case we have modified the VHDL exporting functionality of the gEDA toolsuite to 
generate QHDL. 
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2. Component model specification The QHDL file is tlien parsed to gener- 
ate tlie cireuit expression in wlrielr referenced components appear as symbols. 
This expression is stored in a library file along with information about model 
parameters and the component names of the referenced subcomponents. Note 
that a library file can be treated as a standalone entity for future circuit de- 
signs. When this file is imported at runtime, the referenced subcomponent 
models are dynamically loaded from their respective library files. Now, the 
full (S, L, iJ) parameter triplet can be generated by explicitly evaluating the 
circuit algebra operations. By means of the symbolic operator algebra, the 
final operator matrices and the Hamiltonian are still in fully symbolic form, 
which can be used to generate the quantum master equation or an appro- 
priate stochastic differential equation in symbolic form. This allows for the 
application of analytical model reduction techniques before turning to purely 
numerical methods. 

3. Numerical simulation Define all scalar model parameters and (truncated) 
Hilbert space dimensions, and compute the behavior of the circuit. 

In Table [T] we list the necessary software tools to implement the QHDL circuit 
design workflow. We plan to publicly release our custom tools in the near future. 

2. An example of the QHDL workflow 

In this section, we present a detailed example applying the QHDL workflow to 
the design, analysis, and simulation of an all-optical S'i?-latch as recently proposed 
in [^Q- The elementary component models {{Sj,lij, Hj)} required for this circuit 
are the following: 

1. Beamsplitters 

y^y^smt' cos^ J J 

2. Phase-delays ee (e*'^,0,0) 

3. Coherent displacements VF(a) = {l,a,0) which models a laser source 
outputting a coherent field with amplitude a € (D, 

4. Kerr- nonlinear cavity (here a unidirectional ring cavity with two input / out- 
put ports) 




As the circuits we discuss here are meant to be used as logical gates in larger circuits, 
we need not include the laser sources in our circuit schematics. Instead, they can 
easily be added at the level of the circuit algebra by feeding a concatenated block of 
laser displacements (sources) into the full network Qwith input = Q < [Wai ffl Wa2 H 
• • • ffl W^„). 
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Table 1. List of software components necessary to realize our QHDL-workflow. 



R,equirement 


Our solution 


Alternatives 


Graphical schematic capture 
tool with VHDL/QHDL export 
capabilities 


gsclicm and gnetlist from the 
gEDA suite 


Graphical design tool from sys- 
tem modeling environments / 
modeling languages, such as 
Modclica O 


QHDL-Parser that computes 
the circuit expression 


A custom parser written in 
Python using the open source 
PLY [3 package 


A parser for computing a cir- 
cuit expression from the Mod- 
elica specification, written in, 
e.g. Mathematica 


Symbolic computer algebra 
system with support for: non- 
commutative operator algebra, 
commutative scalar algebra co- 
efficient, operator-valued ma- 
trix algebra and the Gough- 
James circuit algebra 


A custom computer algebra 
system written in Python llSH 
and interfacing with SymPy 
[19|] for the scalar coefficient al- 
gebra 


Mathematica I20II plus an im- 
plementation of the Gough- 
James circuit algebra 


Numerical backend to con- 
vert symbolic operator expres- 
sions into matrices and simu- 
late the system 


Gustom algorithms for solving 
the Master equation as well 
as quantum stochastic differen- 
tial equations implemented in 
Python and G using optimized 
numerical libraries for linear al- 
gebra [2l[2l. 


The Quantum Optics Toolbox 
VM for MATLAB 0| or simi- 
lar library for modeling the dy- 
namics of open quantum sys- 
tems, such as QuTIP 



(a) The two-cavity pseudo-NAND-latch 

We have recently proposed [Ij several different optical circuits to realize three 
classical logic gates: an AND-gate, a NOT-gate with integrated fanout of two and 
a combined (but imperfect) NAND-gate (Figure [3ja)), which in the following we 
will call pseudo-NAND gate as it only works properly when at any given time at 
least one input is in the 'on' state. The first two of these gates used in sequence also 
realizes a NAND gate, but the advantage of the pseudo-NAND is that it requires 
only a single Kerr-nonlinear cavity component. The QHDL workflow can be readily 




Figure 3. Pseudo-NAND circuit schematic (a) as created with gschem and its device 
symbol embedded as a component in a SR-NAND-latch circuit (b). 



applied to design the pseudo-NAND circuit and automatically generate the circuit 
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expression in terms of its componentflj: 



(ll ffl ((li ffl (($ E li) < B2)) < P(i32) < (i^ B li))) 

<(BiE(P(2i)<(VKEli))) 



(2.1) 



where the beamsplitter symbols are defined by Bi = Qbs{j), B2 = Qbs{(^), the 
output correction phase is $ = U{(j)), the constant coherent displacement compo- 
nent is W = W{I3) and the Kerr cavity is given hy K = Qk{^,X^ ki = K2 = k). 
The network expression (|2.ip looks complicated but can be verified easily by com- 
paring its visual representation (Fig. 4(a) with the original circuit schematic. 



Moreover, since the scattering matrix of K is in block-diagonal form, it is possible 
to decompose the cavity component K = K1SK2, where the Hamiltonian of K can 
be assigned to either of the two blocks. Upon substituting this decomposable form 
into the expression, the automatic expression simplification built into our circuit 
algebra implementation yields the following form: 

{ (ll B Ki) < Bi} E { ($ B ll) < ^2 < (W/' B K2) }, (2.2) 

which is visually represented in Figure |4(b)[ The numerical model parameters as 











B, 






A" 




H ... h 




B2 i 















(a) Circuit expression 112.111 as generated by the QHDL-Parser 



(b) Simplified expression 112. 2t 

Figure 4. Pseudo-NAND circuit expression visualizations. As can easily be verified visually, 
the simplified expression follows from decomposing K — K\ ffl K2 and 'pulling' K2 down 
into the fourth row. These expression simpUfications are automatically performed by our 
symbolic circuit algebra software. 



given in are 6 = 0.891, x ~ ~5/6, A = 50, k, = 25, (j) = 2.546 and the auxiliary 
constant input amplitude is given by /3 = —34.289 — 11.909i. The coherent input 
amplitudes corresponding to the logical signals 'on' and 'off' are then given by 
a = 22.6274 and 0, respectively. 

f 1 2 ... n \ . , . . 

I Jdere we represent a permutation cr = I (n)) ^^^^ image tuple 

(o-(l)a(2) ...a{n)). 

t These visualizations were automatically generated using another software tool we have im- 
plemented. 
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Listing 5. QHDL source for the pseudo-NAND latch. 



pseudo—NAND latch 

ENTITY NS3fR3fANDXATCH IS 

PORT (NS, Wl, kerr2_extra, NR, W2, kerrl.extra : in fieldmode; 

BSl.l.out, kerrl_out2, OUT2_2, BS1.2_out, kerr2_out2, OUT2_l : out fieldmode); 
END NS.NR3fANDXATCH; 

ARCHITECTURE nctlist OF NSJVR.NANDXATCH IS 
COMPONENT nand 

PORT (A, B, WJn, kerr_in2 : in fieldmode; 

uol, kerr_outl, NAND_AB, OUT2 : out fieldmode); 
END COMPONENT; 

SIGNAL FB12, FB21 : fieldmode; feedback signals 

BEGIN 

NAND2 : nand 
PORT MAP ( 

A => NR, B => FB12, WJn => W2, kerr_in2 => kerr2_extra, 

uol => BS1.2_out, kerr_outl => kerr2_out2, NAND_AB => FB21, OUT2 => OUT2_2); 

NANDl : nand 
PORT MAP ( 

A => NS, B => FB21, WJn => Wl, kerrjn2 => kerrl.extra, 

uol => BS1.2_out, kerr_outl => kerrl_out2, NAND_AB => FB12, OUT2 => OUT2_l); 
END netlist; 



As in classical circuit theory, two NAND-gatcs in a mutual feedback config- 
uration, as shown in Figure [3jb), can be used to realize a latch with inverted 
inputs S and R. A latch features controllable bistable behavior and thus realizes 
a single-bit memory unit. It has two inputs: S(ET) and R(ESET), which can be 
activated individually to control the internal logical state to 'on' or 'off', respec- 
tively. Ideally, when both S and R are 'off' (HOLD-condition), the internal state 
remains stable. In practice, quantum fluctuations and noisy inputs lead to spon- 
taneous switching between the two internal states. One of the design goals is thus 
to decrease the rate at which this spontaneous switching occurs. The QHDL code 
as produced by gnetlist [l^ (slightly edited to be more concise) can be found in 
Listing [5] and the circuit component library file generated by the QHDL-parser 
is presented in Listing \6\ in [Appendix B[ Substituting the individual component 
models into the circuit expression yields the full triplet (So,Lo,iJo) for the latch. 
Finally, after feeding in the coherent input signals S and R into their respective 
ports: (S, L, i7) = (So, Lq, Hq) < {W{S) E I2 ffl W(R) ffl I2) the parameters assume 
the following form 




(2.3) 
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■ sin 



sin 9e 



9e"f' b - 

KCOSI 

a ^ 



■■ COSf 



_s ^ 

b + 13 sin 6 
,^__^cos0e^^ 



2 ""^"^ ' V2 + x/2 

cos 6 a + P sin 6* 



COSf 



,40 



J4> 



(2.4) 



H = A (a^a + b^b) + x (a^'a^"aa + b^b^bb) -j= sin 6* sin (/) (afe^ + a}l 

v2 



V2k 
H :— « 



S +^*cos( 



h.c. 



/2k 



R +/3*cos6'e" 



(2.5) 
fe- h.c. 



Due to the symmetry of the midcrlying circuit model, the model parameters are 
invariant under exchange of the two pseudo-NAND gates, which corresponds to 
simultaneously exchanging S O R, (a, a^) o (6, 6^), (ii, L2, -^3) ^ {L4, L^^Lq) and 
Si -;-> S2. This symmetry suggests that the most likely candidates for the internal 
logical states 'on' and 'off' correspond to the case where one internal cavity mode 
is in a high power state and the other one in a low power state and the opposite 
case, obtained by exchanging the cavities states. This is indeed the case, and in fact 
it follows from the basic way in which we have designed our pseudo-NAND gate; 
'on' <^ {NANDl cavity power is low, NAND2 cavity power is high} and 'off' <^ 
(NANDl cavity power is high, NAND2 cavity power is low). 

To understand our model's dynamic behavior we turn to numerical methods. 
The simulation of this model is carried out by representing the operators as nu- 
merical matrices in a truncated product basis of Fock-states of total dimension 
= 75^ = 562qj]. We carried out a large number of quantum jump trajectory 
simulations 23|, |26[ with the following sequence of alternating input conditions: 
0.5 time units of SET, 5 units of HOLD, 0.5 units of RESET, 5 units of HOLD 
(repeated twice) . The upper plot in Figure [5] presents a typical simulated trace 
where the system is subjected to this sequence of input conditions. We generally 
find that the SET and RESET input conditions successfully drive the system into 
the desired cavity states, while the cavities remain in their states during the HOLD 
condition. Although a simulation of the full master equation is feasible using current 
HFC hardware and sparse matrix storage [8|, quantum jump simulations exhibit 
the inherently bistable nature of our synthesized latch more clearly. 



(b) Model reduction in the (S,L,ff) context 

As the latch could readily be used as a coniponent in more complex circuits, 
such as flip-flops or even quantum memories 0, @| , it would be highly desirable to 
reduce the Hilbert space dimension A'^^ required to represent it. Since we are work- 
ing with quantum circuit models, we are ultimately limited by exponential scaling 
of the state space with the number of components (although it may be possible to 
develop efficient simulation procedures when components are only weakly entan- 
gled, as should be the case in ultra- low power classical signal processing). However, 

f I.e., each individual cavity basis is given by {|0) , \1) , . . . \N — 1)}. 
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full quantum model 




reduced quantum model 




time 

Figure 5. Simulated input sequence for the full pseudo-NAND latch model (upper) and 
our reduced model (lower). The red trace (lighter, in grayscale) is given by (a^a) and the 
blue trace presents {b^b). The SET and RESET input conditions, marked by the short 
intervals between dashed vertical lines, induce transitions to their respective target latch 
states: 'on' corresponds to the a-mode being in a high photon number state, while 'off' 
corresponds to a high photon number in mode b. 

there is clearly much to be gained by developing accurate model reduction proce- 
dures that allow us to replace high-dimensional ab initio models for components 
within a circuit by much lower-dimensional effective models. Such model-reduction 
strategies could presumably be applied hierarchically. As in the classical theory of 
signals-and-systems there are many potential strategies for dimensional reduction 
of quantum inp ut-output models, although little work has yet been done on this 
subject [23, dg. Here we describe an empirical approach, similar to the classical 
strategy of approximating Markov chains [29^ , which utilizes numerical simulation 
and statistical analysis to derive an effective (S, L, H) model for the pseudo-NAND 
latch suitable for embedding within a more complex circuit. 

Our approach is based on the assumption that the device state can be inferred 
with reasonable accuracy from a small number of observables. We can then con- 
struct a dynamic model just in terms of these parameters [l^l- By generating many 
quantum jump trajectories for the full (S,L,i7) model, we generated many time 
series for the expectation values of the cavity field photon numbers (a^a) and 
under the three valid input conditions {HOLD, SET, RESET}. We can now coarse- 
grain the 2-dimensional space of expectation values and associate an internal model 
state i € {1, 2, . . . , M} with each bin that is actually visited during the trajectory 
simulations, but due to the high correlations between the cavity photon numbers, 
we are actually able to obtain good results by performing this coarse-graining or 
'binning' procedure in terms of the single quantity 

D = (a^a) - (b^b) , 
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implying that within the two-dimensional configuration space our system always 
stays fairly close to a one-dimcnsional submanifold. 

By analyzing the observed transitions between these reduced states for each 
input condition ^ G {HOLD, SET, RESET}, we calculate an empirical estimate 
p(4) = (j>[f)i fj^i of the conditional transition probabilities = P{xn+i = j\xn = 
i, ^) and thus model the system in terms of a discrete time Markov chain with a set 
of conditional transition probabilities for each particular input condition ^. The time 
step 6t of the discrete Markov chain corresponds to the interval at which we sampled 
our original continuous-time system. We now wish to get back to a description that 
is compatible with our (S,L,iJ) formalism. In the following we briefly outline a 
procedure to do this: For a temporally homogeneous Markov jump process with an 
even number of states is {1,2,..., M} and transition rate matri^lj Q = (7ij),f 
we can define a X-channcl (So, Loi -^o) model with states corresponding directly to 



the Markov process states {|1) , |2) , . . . , \M)} via 

So = Ik, (2.6) 

Lo = (V7»iji (ill , ■■■,V^iKjK Uk) {iK\f (2.7) 
Ho = 0, (2.8) 

where the components of Lo drive transitions {ik jk,k = 1,2... K} and K is 



given by the number of positive transition rates jij > 0. By construction, as one 
may verify by writing down the master equation, this system always collapses into 
a purely classical mixture of the coarse-grained states. Equivalently, in a quantum 
jump trajectory simulation, after the first quantum jump, the state is always given 
by a single such state. In fact, in such a trajectory simulation, this system behaves 
exactly like the original Markov jump process. Neglecting for now that our original 
model has three different input conditions ^ G {HOLD, SET, RESET} and thus 
three different conditional transition matrices P^^', we first discuss how to move 
from the discrete time Markov chain model to a continuous time Markov jump 
process. Rephrasing this question, we can ask the following: is there a Markov 
jump process with conditional transition probability matrix P{t) that 'looks like' 
our Markov chain when stroboscopically probed at fixed time intervals 6t7 If our 
Markov chain has transition matrix P, then we need to determine a generator 
matrix Q such that 

P = P(< = St) EE e^*^ « 1 + (5tQ + 0{St^). (2.9) 
If our sample interval dt is sufficiently small, we may define 

Q^l(P-l) (2.10) 

as an approximation to the conditional rate matrix. We now carry out the procedure 
outlined above to create a model (So,Lo,-ffo) that realizes the HOLD condition. 
The transition rates of the HOLD condition alone lead to a system that has two 
bistable clusters of states with low state indices and high state indices, respectively. 

To account for the input-controlled switching in the SET and RESET condi- 
tions, we extend our model by concatenating it with a second model that explicitly 

t Q is also often referred to as the generator matrix. 
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includes the input fields {S^,'Lgj^,H^) = (Si,Li,iJi) <] iW{S) ffl W{R) ffl I2). 
Hence, in the SET and RESET conditions, the HOLD transitions continue, but we 
drive further transitions through this additional component. Here, (Si,Li,_ffi) is 
given by 



Si 

Li 

Hi 



/4ss 








R 


-Si? 







\ 







-a{l 



1 



0, 



(2.11) 

(2.12) 
(2.13) 



where the 'drift' operators E5 and are defined by 

Ss = \M - 4) {M - 1| + |M - 6) {M - 3| + • • • + |M/2 - 1) (A//2 + 2| , 
Sfl = |5) (2| + |7) (4| + • • • + \AI/2 + 2) {M/2 - 1| . 

For our simulation we chose M = 38 <C iV^, but the general ansatz works for 
a range of different AI ^ 4k + 2 with sufficiently large k. The drift operators 



satisfy S^y^ 



= 0, 



as well as the projection relations (E^^^S^/^)^ = E^^^E^/g, 

^R/S^R/s f = ^R/s^R/s and Efl/5E^y_5Efl/5 = E^/g. These relations suffice 
to show that Si as defined above is indeed unitary. To make sense of the effect 
of this extension to our model, consider now what happens for the different input 
conditions: In the HOLD condition S = R = a the input fields cancel out all 
elements of the coupling vector Li and we have (Sgjj, Lgjj, i/gjj) = (Si, 0,0), i.e. 
the transition dynamics of our full system (Sg^;, Lg^;, _Hgp;)ffl(So, Loi ^0) are simply 
given by those of (So, Lq, iJo) alone. 




Figure 6. Here we schematically visualize the state space and the transitions of the reduced 
model. The SET transitions (red, from right to left) introduce a drift that drives the 
system to the states on the left, corresponding to the logical 'ON' state of the latch. 
The RESET transitions (blue, from left to right) have the opposite effect. The HOLD 
transitions (depicted in grey) are always active, but in the absence of additional SET and 
RESET transitions only very rarely lead to a switch of the logical latch state. 

In the SET condition, however, we have S = 0, R = a and thus 

(SsR,Lsj^,i/5R) = (Si,(-a(l-E^<,E5), 0, -aE5,0)^,o). 

The full system (Sg^-, Lgp-, _Hgp-) ffl (So,Lo,iJo) now features the drift operator as 
an additional transition operatoiQ — aEs which induces transitions {M — 1 ^ 

f The second non-zero element of Lgj^, which is a projection operator, does not affect the 
transition dynamics due to the fact that our system is never in a superposition of two states. 
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M-4, M-3 ^- M-6, . . . ,M/2 + 2 ^ M/2-1} with constant rate |ap. Together 
with the HOLD transitions, these lead to a drift from states with high index (cor- 
responding to the logical 'off' state of the latch) to those with low index ('on'). On 
the other hand, in the RESET condition, the situation is reversed. Now the other 
transition operator of Li is canceled out and the non-zero transition operator — aEj^ 
drives transitions in the reverse direction {2 — > 5, 4 — > 7, . . . , AI/2 — 1 — > Af/2 + 2}, 
again with constant rate jap. In Figure IH] we visualize the transition structure of 
the model schematically. 

Note also that we can emulate the state-dependent coherent output field(s) of 
the latch by concatenating a triplet (Sout, Lout, -ffout) that re-routes bias input fields 
via state-dependent scattering into one or more output channels. For example we 
could use (Sgp;, Lgp;, Hg^) ffl (So, Lq, i?o) H (Sout, Lout, 0) where 

I -x /-I /e*"^" cosft; — e*'^"sin0A , ^ ^,7^ /„ -, ,x 

Sout = El^)(^l(,*.3in^. e*-coseJ' ^out = Sout(/3', 0)^, (2.14) 

where /3' is the complex amplitude of a bias field and the parameters {di, (pu, 02^} 
arc chosen such that the outputs of (Sout, Lout, 0) vary as desired with the internal 
state Having thus created a reduced model that mimics the desired input- 
output behavior in (S, L, H) form, we can use it to replace the full latch model in 
more complex circuits. If we had already specified a QHDL file for such a circuit, 
we could simply replac^ the referenced latch component with the reduced model 
component. Re-parsing this modified QHDL-file would then yield a computationally 
more tractable model for simulations. 

3. Conclusion 

In this paper we have described the use of QHDL, a quantum hardware description 
language, to facilitate the analysis, design, and simulation of complex networks con- 
structed from interconnected quantum optical components. We have also presented 
a parsing algorithm for obtaining quantum equations of motion from the QHDL 
description. QHDL can be used as the basis for a schematic capture workflow for 
designing quantum circuits that automates many of the conceptually challenging 
and computationally demanding aspects of quantum network synthesis. As QHDL 
inherits the hierarchical structure of VHDL, its use may facilitate the crucial devel- 
opment of hierarchical model reduction methods for quantum nonlinear photonics. 

Important future directions for QHDL research include simulation strategies 
for exploiting weak entanglement among components, stability analysis and design 
optimization of QHDL-based models [31| , and the incorporation of techniques from 
static program analysis and formal verification to assist in the design of complex, 
hierarchically defined photonic components. While we have emphasized classical 
photonic logic as a tutorial paradigm for QHDL in this paper, emerging ideas in 
quantum information processing and quantum sensing/metrology may provide even 

f In principle it should be possible to include the reduced model as an alternative architecture 
for the latch entity and to select whether or not to use it in place of the full model at compile-time 
using a VHDL configuration file. However this would require some enhancements to the QHDL- 
Parser to correctly handle the K extra (vacuum) input ports required by the reduced model to 
drive spontaneous transitions among the internal states. 
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more compelling applications for QHDL as a convenient and extensible modeling 
framework. 

This work is supported by DARPA-MTO under Award No. N66001-11- 1-4106, by tlie 
National Science Foundation under Grant No. PHY-1005386, and by the SU2P Program 
of Stanford University and Research Councils UK. 



Appendix A. Reduced parameters in case of signal feedback 

Upon feeding the fc-th output channel of a system Q ~ {8,1,, H) back into its i-th input, we 
get a system [ (S, L, H) = ^S, L, fl^ with one less channel cdim [ Q ]^.^; = cdimQ — 1, 
where effective parameters are then given by P| 

/ Si, \ 

S21 



s = s 



Sk-ii 
Sk+1 1 

\ Snl J 



(1 — Ski) ^ {Ski Sk2 



Ski+1 



(Al) 



+ 



S21 



Sk-i I 
Sk+1 1 



\ S„i I 



{\-SkiY^Lk 



H ^ H + 'S 



U=i 



{i-Skiy^Lk 



(A 2) 



Here we have written as a shorthand notation for the matrix S with the fc-th row 

and l-th column removed and similarly L^,^ is the vector L with its k-th entry removed. 
These resulting parameters fulfill the condition^ for circuit components. Moreover, they 
have shown that in the case of multiple feedback loops, the result is independent of the 
order in which the feedback operation is appliecjj. 



Appendix B. Latch circuit library file 

Listing 6. Python[l^ source for the pseudo-NAND latch circuit library component. 
#!/usr/bin/env python 

from qhdl_component_lib. library import retricve_component , ma]ie_namespace_string 
from qnet.qos_algebra import P_sigma, qid, FB 
from sympy. core. symbol import symbols 

t This is obvious for L and H, for a proof that S is indeed unitary see Gough and James's 
original paper 

I Note however that some care has to be taken with the indices of the feedback channels when 
permuting the feedback operation. 
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NCHANNELS = 6 
GENERIC.DEFAULT.VALUES = {} 

def arch_default(name_space = ", **generic_params): 

Generate a symbolic LATCH expression using the provided name^space. 
Return a 2— tuple (netlist^symbolic, var^map) 
such that calling 

>>> SLH = netlist^symbolic.substitute(var^map).evalf() 
will result in an SLH triplet, where the parameters are accessible via 

»> SLH.S 

»> SLH.L 

»> SLH.H 

# dictionary that stores the replacement SLH models for 

# the circuit component symbols 
var_map = {} 

# load symbolic component expressions as well 

# as actual SLH model replacements for NANDl... 
NANDl, NANDl_var_map = retrieve_component('nand', 4, 

make_namespace_string(name_space, 'NANDl')) 
var_map.update(NANDl_var_map) 

#...and NAND2 

NAND2, NAND2_var_map = retrieve_componcnt('nand', 4, 

make_namespacc_string(name_space,'NAND2')) 
var_map.update(NAND2_var_map) 

########## symbolic circuit expression as computed by the parser 

# the '+' operator is overloaded to the concatenation— operation [+] 

# the '<<' operator is overloaded to the series— operation <| 

# P^sigma(s(0),s(l),...(s(n—1))) is a zero— based channel permutation object 

# FB(Q, k,l) is the feedback operation [Q]^{k—>1} (with zero— based channel indices) 

# qid(n) is the n— channel identity 
netlist_symbolic = (P_sigma(0, 2, 3, 4, 5, 1) 

<< FB( 
( 

(qid(l) + 

({qid(3) + P.sigma(2, 0, 1)) 

<< ((P^igma(l, 2, 3, 0) << NAND2) + qid(2)))) 
<< P_sigma(0, 2, 5, 6, 1, 3, 4) 
<< ((P^igma(0, 3, 1, 2) « NANDl) + qid(3)) 
), 6, 1) 

« Pjigma(0, 4, 5, 3, 1, 2)) 
return nctlist_symbolic, var_map 



References 

[1] J. Gough and M. R. James. Quantum Feedback Networks: Hamiltonian Formulation. 
Communications in Mathematical Physics, 287(3):1109-1132, December 2008. ISSN 
0010-3616. doi: 10.1007/s00220-008-0698-8. 

[2] John Gough and M.R. James. The Series Product and Its Application to Quantum 
Feedforward and Feedback Networks. Automatic Control, IEEE Transactions on, 54 
(ll):2530-2544, June 2009. ISSN 0018-9286. doi: 10. 1109/TAC. 2009. 2031205. 



Article submitted to Royal Society 



Specification of photonic circuits using QHDL 



19 



[3] H.J. Carmichael. Quantum Trajectory Theory for Cascaded Open Systems. Physical 
Review Letters, 70(15):2273-2276, Aug 1993. doi: 10.1103/PhysRevLett. 70.2273. 

[4] CW Gardiner. Driving a quantum system with the output field from another 
driven quantum system. Physical review letters, 70(15):2269-2272, 1993. doi: 
10.1103/PhysRevLett.70.2269. 

[5] Joseph KerckhofT, Hendra Nurdin, Dmitri Pavhchin, and Hideo Mabuchi. Design- 
ing quantum memories with embedded controh Photonic circuits for autonomous 
quantum error correction. Physical Review Letters, 105(4) :040502, jul 2010. ISSN 
0031-9007. doi: 10.1103/PhysRevLett. 105. 040502. 

[6] J Kerckhoff, D S Pavhchin, H Chalabi, and H Mabuchi. Design of nanophotonic 
circuits for autonomous subsystem quantum error correction. New Journal of Physics, 
13(5):055022, May 2011. ISSN 1367-2630. doi: 10.1088/1367-2630/13/5/055022. 

[7] Hideo Mabuchi. Coherent-feedback control strategy to suppress spontaneous switch- 
ing in ultralow power optical bistability. Applied Physics Letters, 98(19):193109, 2011. 
ISSN 00036951. doi: 10.1063/1.3589994. 

[8] Hideo Mabuchi. Nonlinear interferometry approach to photonic sequential logic. Ap- 
plied Physics Letters, 99:153103, 2011. 

[9] Kurt Jacobs, H.I. Nurdin, F.W. Strauch, and Matthew James. Frequency Conversion: 
Side-band cooling, state-swapping, and coherent control of mechanical resonators. 
Arxiv preprint arXiv:1003.2653, pages 3-6, 2010. 

[10] J. E. Gough, M. R. James, and H. I. Nurdin. Squeezing components in linear quantum 
feedback networks. Physical Review A, 81(2), February 2010. ISSN 1050-2947. doi: 
10. 1103/PhysRevA. 81. 023804. 

[11] C.W. Gardiner and P. ZoUer. Quantun Noise: A Handbook of Markovian and 
Non-Markovian Quantum Stochastic Methods with Applications to Quantum Optics. 
Springer Series in Synergetics, 2nd edition, 2000. 

[12] L. Bouten, R. Van Handel, and M. James. An introduction to quantum filtering. 
SIAM Journal on Control and Optimization, 46(6):2199-2241, 2007. doi: 10.1137/ 
060651239. 

[13] Howard M. Wiseman and Gerard J. Milburn. Quantum Measurement and Control. 
Cambridge University Press, Cambridge, 1st edition, 2010. 

[14] Volnei A. Pedroni. Circuit Design and Simulation with VHDL. MIT Press, Cam- 
bridge, 2004. 

[15] The gEDA project. gEDA - gpl electronic design automation suite, 1998 -. URL 
|http : //www ■ gpleda . org[ 

[16] The Modelica Association. Modelica, 1996-. URL [http : / /modelica . org[ 

[17] David Beazley. Ply - python lex-yacc, 2001 -. URL |http : //www . dabeaz . com/ply] 

[18] Guido Van Rossum. The python language, 1990-. URL [http : //www . python . org/ [ 

[19] SymPy Development Team. Sympy - python library for symbolic mathematics, 2007 
-. URL [http: //www, sympy ■ org [ 



Article submitted to Royal Society 



20 N. Tezak, A. Niederberger, D. S. Pavlichin, G. Sarma, and H. Mabuchi 



[20] Wolfram Research, Inc. Mathematica, 1988 -. URL 

: / / www .wolfram.com/mathematica 

[21] David Ascher, Paul F. Dubois, Konrad Hinsen, James Hugunin, Eric Jones, Travis 
Oliphant, and Pearu Peterson. SciPy and Numpy: Open source scientific tools for 
Python, 1996-. URL http: //www. scipy . org/, 

[22] Intel Corporation. MKL - Math Kernel Library, 2003-. URL 

: //sof twajre . Intel . com/en-us/ intel-nikl| 

[23] Sze M Tan. A computational toolbox for quantum and atomic optics. Journal of 
Optics B: Quantum and Semiclassical Optics, 1(4):424, 1999. 

[24] The MathWorks, Inc. MATLAB - the language of technical computing, 1984 -. URL 
|http : //www.mathworks . com/products/matlab| 

[25] J.R. Johansson, P. D. Nation, and Franco Nori. QuTIP: An open-source Python 
framework for the dynamics of open quantum systems. AnAv preprint arXiv.lllO. 057, 
2011. URL Ettp : / /code . google . com/p/qutip] 

[26] S. Haroche and J. M Raimond. Exploring the quantum : atoms, cavities and photons. 
Oxford University Press, Oxford; New York, 2006. ISBN 9780198509141; 0198509146. 

[27] L Bouten, R van Handel, and A Silberfarb. Approximation and limit theorems for 
quantum stochastic models with unbounded coefficients. Journal of Functional Anal- 
ysis, 254(12):3123-3147, June 2008. ISSN 00221236. doi: 10.1016/j.jfa.2008.02.013. 

[28] J.E. Gough, H.I. Nurdin, and S. Wildfeuer. Commutativity of the adiabatic elimina- 
tion limit of fast oscillatory components and the instantaneous feedback limit in quan- 
tum feedback networks. J. Math. Phys., 51(12):123518, 2010. doi: 10.1063/1.3520513. 

[29] H.J. Kushner and P. Dupuis. Numerical Methods for Stochastic Control Problems in 
Continuous Time. Springer- Verlag, New York, 2001; 1992. 

[30] A.E.B. Nielsen, Hopkins A.S., and Mabuchi H. Quantum filter reduction for 
measurement-feedback control via unsupervised manifold learning. New J. Phys., 
11:105043, 2009. 

[31] Armand Niederberger, Nicolas R. Ganger, and Hideo Mabuchi. Efficient quantum 
design optimization and stability analysis, in preparation, 2011. 



Article submitted to Royal Society 



